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Estimation de R xx 
sur K' echantillons 



Calcul de R: 



ft Precede et dispositif permettant a un modem de se synchroniser sur un transmetteur de donnees 
numeriques par voie hertzienne en presence de brouilleurs. 

(St) Le procede consiste & estimer (9) la matrice 
de correlation ftxx des signaux recus du trans- 
metteur sur un ensemble de n capteurs de 
reception, & calculer (10) la matrice de correla- 
tion inverse ft£, a calculer (11) des vecteurs 
d'intercorrelation t m © ntre ,es signaux recus 
sur Tensemble des n capteurs et un signal de 
replique d(k) connu, & calculer un critere de 
synchronisation (15) en effectuant les produits 
des vecteurs d'intercorrelation transposes 
conjugues ft^ et de la matrice des vecteurs 
d'intercorrelation t m * et comparer (16) la 
valeur du critere obtenu & une valeur de seuil t\ 

determinee pour placer la synchronisation sur 

Techantillon pour lequel la valeur des criteres 

depasse la valeur de seuil y\. Pour estimer le 

decalage Doppler, la sequence de synchronisa- 
tion est* decoupee en tranches d'un nombre 

determine de symboles pour effectuer des cal- 

culs de transformation de Fourier rapide sur 

des vecteurs d'intercorrelation determines sur 

chaque tranche. 
Application : Transmission radio de donnees. 
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les vecteurs r Kdi 



Calcul de fFT^ sur 
les vecteurs r x * y de 
tales T xd ^4 ' 



If) 

CO 
CD 

m 

CO 
CO 

o 




estimation decalag 
en frequence 
|lca!cufde TFD) 



FIG. 7 



EP 0 665 665 A1 



' La presente invention cencerne un procede et un dispositif perrnettant a un modem de se synchroniser 
sur un transmetteur de donnees numeriques par voie hertzienne notamment en presence de brouilleurs. 

Dans les systemes de transmissions de donnees numeriques par voie hertzienne la synchronisation est 
effectuee, k partir de sequences de symboles connus ou recepteur et intercalees entre les symboles d'infor- 
5 mation, par correlation entre le signal recu et un signal de reference construit a partir des symboles connus. Le 
resultat est compare k un seuil pour decider de la synchronisation. Des systemes perrnettant la synchronisation 
d'un modem k partir d'un seul capteur et d'un seul recepteur sont par exempie connus de I'article D.V. Sarwate et 
M.B. Pursely intitule "Cross correlation properties of pseudo random and related sequences" paru dans IEEE Vol. 
66, pp. 593-619, May 1980. Mais les techniques mises en oeuvre qui perrnettent en outre une estimation des,ca- 
10 racteristiques du canal de transmission deviennent rapidement inef f icaces en presence de brouillage. 

Lorsque le rapport signal k bruit n'est pas suffisant pour la prise de synchronisation par une technique 
mono-capteur, des techniques visant k accroitrele rapport signal k bruit et a permettre la synchronisation sont 
utiiisees. 

Ces techniques peuvent etre regroupees en deux families. Une premiere famille vise k pretraiter le signal 
15 recu avant d'effectuer la synchronisation du signal pretraite. La seconde famille vise k ameliorer le rapport 
signal k bruit et k effectuer la synchronisation conjointement. 

Dans les techniques relevant de fa premiere famille il faut egalement distinguer deux categories : 

- Une premiere categorie met en oeuvre un pretraitement mono-capteur de type excision de brouilleurs 
avant d'effectuer la synchronisation. Ce traitement est connu par exempie de I'articie de LB. MILSTEIN 

20 "Interference rejection Techniques in spread spectrum communications", paru dans Proc. IEEE, Vol. 76, 

pp. 657-661, June 1988. Suivant cet article on cherche k predire la contribution du brouilleur, dont la 
bande est supposee plus petite que la bande du signal, sans predire celle du signal utile. Cela permet 
d'ameliorer les performances de la synchronisation en presence d'interferences bande etroite (type 
PSK). Mais le traitement reste inefficace en presence d'un brouillage occupant toute la bande utilisee 

25 par le signal utile. Cela est en partie due au fait que les techniques de pretraitement mono-capteur ne 

perrnettent un traitement des signaux que dans un espace limite a I'axe des temps et k I'axe des fre- 
quences. 

- La seconde categorie concerne les techniques mettant en oeuvre un pretraitement multi-capteurs de 
rejection, precedant la prise de synchronisation. L'idee consiste a mettre en oeuvre un f iltrage spatial, 

30 discriminant le signal et les brouilleurs par leur direction d'arrivee, le temps, le spectre ou la puissance 

comrrie cela est decrit dans I'ouvrage de R.A. Monzingo, T.W. Miller, intitule "Introduction to adaptive 
arrays", edite par John Wiley and Sons, New-York, 1 980 et la these de P. Chevalier, intitulee "Antenne 
adaptative : d'une structure lineaire k une structure non lineaire de Volterra", These de doctorat, Uni- 
versity de Paris Sud, Juin 1991 puis de correier le signal de reference avec la sortie de I'antenne. Les 
35 techniques de f iltrage spatial discriminant le signal et les brouilleurs par lo temps, le spectre ou la puis- 

sance sont generalement sous-optimales, alors que ceiles exploitant la connaissance de la direction 
d'arrivee du signal necessitent la mise en oeuvre d'une technique de goniometrie, couteuse et delicate 
dans la gamme des frequences HF notamment du fait des multi-trajets. 
Pour ces raisons, les techniques de la seconde famille, visant a effectuer la rejection de brouilleurs conjoin- 
40 tement avec la prise de synchronisation, a partir de la connaissance du signal de reference, apparaissent plus 
prometteuses. Une telle technique a par exempie ete proposee par M. Compton dans son article "An adaptive 
array in a spread spectrum communication system" paru dans Proc. IEEE, vol. 66, pp. 289-298, March 1978, 
pour la synchronisation de signaux PN, k partir d'une architecture analogique integrant des "Delay Lock Loop". 
Une seconde technique fut proposee par L.E. Brennan et I.S. Reed dans leur article intitule "An adaptive 
45 " array signal processing algorithm for communications", paru dans IEEE Trans-Aero-Syst. vol. AES-18, n°1, 
pp. 1 24-1 30, June 1982 pour les transmissions numeriques et completee par D.M. Duglos et R.A. Scholtz dans 
I'article intitule "Acquisition of spread spectrum signals by an adaptive array", paru dans IEEE Trans-Acous- 
Speech-Signal-Proc, vol. ASSP-37, n°8. pp. 1253-1270. August 1989. Ces techniques apparaissent teutefois 
beaucoup plus couteuses que les techniques mono-capteur et sont malheureusement inapplicables telles quel- 
50 les en presence des derives Doppler du signal inherentes aux transmissions ionospheriques. 
Le but de I'invention est de pallier les inconvenients precites. 

A cet effet, I'invention a pour objet un procede perrnettant k un modem de se synchroniser sur un trans- 
metteur de signaux numeriques par voie hertzienne en presence de brouilleurs consistant k: 

- estimer la matrice de correlation Rxx de signaux recus sur un ensemble de N capteurs de reception par 
55 blocs composes d'un nombre d'echantilions d'une sequence de synchronisation transmise par le trans- 
metteur, A 

- k calculer la matrice de correlation inverse R»\ 



2 



I 
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- k calculer des vecteurs d'intercorr6lation f Xd entre les signaux X(k) recus sur I'ensemble des N capteurs 
et un signal de replique d(k) connu, 

- k calculer un entire de synchronisation effectuant le produit scalaire entre le yecteur d'intercorreiation 
rxd et le vecteur obtenu par produit de I'inverse de la matrice de correlation R» et du vecteur d'inter- 

5 correlation rxd 

- k comparer la valeur du critere obtenu k une valeur de seuil ti determinee pour placer la synchronisation 
sur rechantillon pour lequel la valeur des criteres depasse la valeur de seuil n, caract6ris6 en ce qu'il 
consiste pour estimer le decalage Doppler inherent aux communications entre 6metteurs et r6cepteurs, 
k decouper la sequence de synchronisation en un.nombre determine de M tranches de ko symboles, k 

10 effectuer sur ces tranches des calculs de vecteurs d'intercorreiation partielles f Xd j (i = 1 M) de telle 

sorte que sur chacune des tranches de ko symboles, le decalage Doppler fasse tourner la phase d'une 
valeur interieure k 90°, k effectuer des calculs de transform^ de Fourier rapide (FFT) sur les vecteurs 
?xdi obtenus, et k calculer une valeur du critere de synchronisation k partir des vecteurs r Xd (kAf), obtenus 
en rdsultatdes calculs de transform^ de Fourier rapide, et de la matrice de correlation Rxx- 
15 L'invention a 6galement pour objet un dispositif pour la mise en oeuvre du proc6d6. 

Le proc6de et le dispositif selon ('invention ont pour avantages qu'ils permettent k la fois d'assurer la syn- 
chronisation d'un modem dans des conditions debrouillage s6v6re.et en presence de multi-trajets, d'estimer 
la derive en frequence entre emetteur et r6cepteur et de detecterles instants d*arriv£e de tous les multi-trajets 
assoctes au signal utile, pourvu que ceux-ci sdient de puissance suff isamment forte par rapport k la puissance 
20 de bruit de fond. 

D'autres caracteristiques et avantages de I'invention apparattront dans la description qui suit faite au re- 
gard des figures des dessins annexes qui represented : 

- La figure 1 un schema fonctionnel d'une antenne adaptative selon invention. 

- La figure 2 un exemple de mise en oeuvre d'un f iltre adapte spatial k discrimination par la modulation 
25 pour la realisation d'une antenne adaptative selon I'invention. 

- La figure 3 un schema d'antenne adaptative pour la mise en oeuvre du procede de synchronisation spa- 
tiale selon I'invention. 

- La figure 4 un graphe pour illustrer la methode d'estimation de la matrice de correlation de signaux pro- 
venant de N capteurs. 

30 - La figure 5 un graphe pour illustrer le d6coupage des symboles de la sequence de synchronisation par 
tranche pour le calcul de vecteur d'intercor relation en presence de decalage Doppler. 

- La figure 6 un graphe montrant revolution de critere de synchronisation mis en oeuvre par I'invention 
en presence d'un decalage Doppler. 

- La figure 7 un organigramme retracant I'ensemble des etapes n6cessaires k la mise en oeuvre du pro- 
35 c£d6 de synchronisation selon I'invention. 

Le procede de synchronisation mis en oeuvre par I'invention fait appel aux techniques connues de syn- 
chronisation mono-capteur pour les appliquer k la realisation d'aritennes adaptatives sur r6plique. 

Selon la technique de synchronisation mono-capteur, le critere de synchronisation C(n) est obtenu en cal- 
culant la correlation y(n) entre K symboles connus dans une sequence de synchronisation et les K echantillons 
' 40 quisuiventl'arriveed'un6chantillonn surlecapteur,en normalisantles resultatsparla puissance p(n)du signal 
recu. Ces calculs sont effectu6s en appliquant les relations : 

45 y(n) = ±£x(k + n)cr(k) (1) 

K k=l 

1 K 

p(n) = ^x(k + -n)x*(k + n) (2) 

55 ou: 

- x(k) est le signal issu du capteur utilise, 

- d(k) est le signal de r6plique connu (K symboles places en debut de trame), 
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- * represents I'operation de conjugaison sur les nombres complexes. 

Ces calculs sont r6iter6s pour chaque echantillon n, et lorsque le critfcre C(n) depasse un certain seuil x\ t 
dependant de la probability de fausse alarme, la synchronisation est decidee. L'instant de synchronisation est 
determine par le rang d'echantillon n 0 pour lequel C(no)> n- La normalisation par la puissance p(n) permet de 
5 choisir un seuil independant du niveau de signal recu. 

Dans ce procede la synchronisation est obtenue lorsque le rapport signal a bruit est suff isant Ses per- 
formances dependent en partie du nombre d'echantillons pris en compte pour ia correlation. Cependant en 
presence de brouilleurs forts, la synchronisation ne peutetre effectuee. 

L*id6e mise en oeuvre dans I'invention consiste alors a effectuer la synchronisation, en sortie d'un traite- 
10 ment d'antenne adaptative, permettant la rejection des brouilleurs. 

Le schema fonctionnel d'une antenne adaptative est represente a la figure 1 , celle-ci comporte un reseau 
1 de N capteurs 1 t ... 1 N , un reseau de filtres adaptes 2 (constitu6s par des circuits de ponderations amplitude- 
phase, desf iltres& r6ponse impulsionnelle f inie RIF ou infinie Rll) couple a un processeur de commande adap- 
tatif 3. Le processeur 3 calcule les coefficients Wi a W N du reseau de filtres 2 pour I'adapter au signal utile 
15 et aux brouilleurs, en I'absence de toute connaissance a priori sur les brouilleurs et a partir d'une information 
minimale sur le signal utile, Les coefficients a W N permettent de former un trou dans le diagramme de 
rayonnement du r6seau d'antenne dans la direction des brouilleurs, tout en maintenant un certain gain dans 
la direction du signal utile. 

Diff6rents types d'information sur le signal utile sont exploites par I'antenne adaptative. Ceux-ci sont d'or- 
20 dre spatial lorsque la direction d'arrivee du signal utile est connue, d'ordre temporelle lorsque les instants de 
presence et d'absence de signal utile sont connus par exemple en evasion de frequence, relatifs a la forme 
d'onde lorsqu'une information sur la modulation du signal utile estdisponible pour disposer d'une replique du 
signal utile en intercalant par exemple des sequences connues dans la trame utile, d'ordre spectral lorsque 
le signal utile occupe une certaine bande de frequence et que les brouilleurs occupentune bande plus impor- 
25 tante, ou statistiques pour s6parer les diff6rentes sources arrivant sur le r6seau en exploitant par exemple les 
statistiques d'ordre 4 des signaux en supposant que les differents signaux sontstatistiquement independants. 

Dans le cas de I'antenne adaptative representee a la figure 2, 1'information sur le signal utile est constituee 
par les sequences de symboles connus intercaiees dans la trame entre les sequences de symboles d'infor- 
mation. Les filtres places en sortie de chacun des capteurs 1 1 a 1 N sont formes respectivement par des circuits 
30 de ponderation amplitude-phase W, a W N des signaux x^tjax^tyfournis par les capteurs. Les signaux obtenus 
en sortie des circuits de pond6ration sont appliques a un circuit de sommation 4. 

Les caracteristiques du reseau de f iltreS 2 sont calcuiees en minimisant une ErreurQuadratique Moyenne 
(EQM) estim6e £ entre un signal de r6plique d(k) et le signal y(k) sortantde I'antenne de la facon sch6matisee 
& la figure 2 oil les elements homologues a ceux de la figure 1 sont reper6s avec les memes references, \ est 
35 def inie par la relation 



40 



4=Z|y(k)-d{k)| 2 =2 w*x(k)-d(kf (4) 
k=l k=l 1 



ou: 



- y(k) est la sortie de I'antenne adaptative et est def inie par 

y(k) = WX(k) (5) 

45 - W est le vecteur poids def inissant I'antenne adaptative, W etant le vecteur transpose conjugue corres- 

pondent, (le signe t en exposant def init I'operateur correspondent.) 

- X(k) est le vecteur forme par les signaux recus sur les capteurs. 

- d(k) est le signal de replique connu. 

La difference y(k) - d(k) est obtenue sur la figure 2 a la sortie d'un circuit soustracteur 5. 
so Le vecteur poids W est calcuie par la formule connue de Wiener : 

W = R-r Xd (6) 

avec : 



55 A 1 ^ t 



Rxx-ijlxOWOO (7) 



4 



EP 0 665 665 A1 

matrice de correlation des signaux provenant des capteurs 1 , a 1 N . (X'(k) etant le vecteur transpose conjugue 
du vecteurX(k)). 



K k=1 



vecteur d'intercorrelation entre X(t) et d(t) (8). 

L'antenne adaptative ainsi utilisee conduit a une maximisation du rapport signal/(bruit + interference) 
connu sous la designation SNIR abreviation anglo-saxonne de "Signal to Noise + Interference Ratio" obtenu 
en sortie de l'antenne, avec une possibility de rejection pouvant aller jusqu'a N-1 brouilleurs, et est appelee 
Filtre Adapte Spatial (FAS). 

La recherche de la synchronisation est effectuee selon I'invention en sortie de l'antenne adaptative en 
discriminant le signal selon sa modulation, et en utilisant le critere de synchronisation mono-capteur. Le sche- 
ma de principe de la synchronisation spatiale est represents a la figure 3 ou les elements homologues a ceux 
des figures 1 et 2 sont represents avec les memes references. Le critere de synchronisation devient en uti- 
lisant les notations precedentes : 



C(n) = l|y(k + n)d(k) (9) 
K k=l 

y(k+nHepresente la sortie de l'antenne adaptative conduisant a une minimisation de I'erreur quadratique 
moyenne \(n) estimee sur K echantillons entre la replique d(k) et le signal arrivant pris a partir de I'echantillon 
n, suivant la relation 

£(n)= ]T|y(k+n)-d(k)| 2 avec y(k + n) = W(n) l .X(k + n) (10) ; 
k=1 

avec 

W(n) = Rxx" 1 (n).fxd(n) (11) 
Les differentes correlations sont calculees sur les K echantillons qui suivent I'arrivee d'un echantillon n 
suivant les relations : 

Rxx(n) = ^ix(k + n)X t (k + n) (12) 
K k=l 



>kd(n)=7ix(k + n)d"(k) (13) 



Lorsque I'echantillon n correspond a la position de synchronisation, les poids W(n) definissent le filtre 
adapte, et y(k+n) contient alors le signal utile d6barrasse de la contribution des brouilleurs. Dans le cas contrai- 
re, aucun signal n'est correle avec la replique d(k) et y(k+n) ne contient plus alors que la contribution du bruit 
de fond. 

L'expression du critere de synchronisation spatiale se simplifie en ecrivant : 



K ( K "\ 

C(n) = -1 2W(n) t X(k+n)d # (k) = W(n) t ± £x(k + n) d*(k) 
K k=1 l K k=1 J 



(14) 
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soit A 

C(n) = ^(n)Rxx(n)-if Xd (n) (15) 

De m§me que pourle cas mono-capteur, ce calcul estr6iter6 pour chaque 6chantillon n, et lorsque le entire 
C(n) depasse un certain seuil n. dependant de la probability de fausse alarme, la synchronisation est decid6e 
(ia normalisation par la puissance est ici obtenue par I'intermediaire de R«). 

L'instant de synchronisation est alors determine par le rang de rechantillon pour lequel C(n 0 ) est superieur 
a une valeur de seuil d6terminee n.- 

Les calculs correspondant sont effectues sur la figure 3 par le bloc 7 de calcul du critere C(n) et par le 
d6tecteur de seuil 8, & partir de la sequence de synchronisation d(t) et du signal y(t) correies par le correlates 
6. 

Cela necessite pour chaque echantillon, d'une part de calculer N correlations entre le signal de r6plique 
d(k) et les signaux x(k) regus par les N capteurs (soit N fois les calculs ef fectu6s pour la synchronisation mono- 
capteur), et d'autre part d'estimer puis d'inverser la rnatrice d'intercorr6lation des signaux capteurs. 

Au cours de ces calculs la synchronisation et le calcul du filtre adapte 2 sont effectu6s conjointement. 
15 Lorsque la valeur du critere depasse le seuil ti, calcuie. la synchronisation est d6cid6e et le vecteur poids de- 
finissant le filtre 2 s'en deduit facilement par la formule connue de Wiener. Le vecteur poids W peut ensuite 
etre utilise pour f iltrer les echantillons contenant Information. 

En supposant que la replique d(t) est de puissance 1 , les correlations tendent vers les valeurs suivantes, 
lorsque le nombre d'echantillons K devient tr£s grand ; 



10 



20 



A 

Rxx — R + 5 S* 



ou: 



- R est la rnatrice de correlation des interferences et du bruit de fond, 
25 - S est la rnatrice du vecteur directeur du signal utile, 

- K 8 est la puissance du signal utile repu sur un capteur. 

Le critere de synchronisation qui converge alors vers une valeur C 0 s'6crit encore: 

_ «. S*R-* S _ SNIRopt m) 
0 " 1 + «. S*R- 1 S " 1 + SNIRopt 1 ' 
30 ou SNIRopt est le SNiR asymptotique obtenu en sortie du filtre spatial adapte 2. 

Lorsque le SNIR est fort, le critere tend vet's 1 , lorsqu'il est faible, il tend vers 0. L'importance du choix du 
seuil tj, fixe la probability de detection de la synchronisation. Le seuil depend naturellementde la probabilite 
de fausse alarme souhaitee. 

Le calcul du critere de synchronisation C, de f Xd et de Rxx est effectue pour chaque position de la syn- 
35 chronisation et notamment la rnatrice Rxx est estimee et inversee egalement pour chaque position de la syn- 
chronisation, la puissance de calcul demandee est £lev6e. Mais il est possible de rgduire la puissance de calcul 
en estimant la rnatrice de correlation Rxx moins souvent sur un intervalle de temps plus important (K* echan- 
tillons au lieu de K avec K'>K de la fagon representee £ la figure 4) et en utilisant cette estimation de Rxx sur 
toutes les positions de synchronisation incluses dans cet intervalle de temps. 
40 Dans ces conditions la rnatrice Rxx s'6crit : 



45 



R XX = 



K' 



K' K' 
2B(t k )B t (t k ) + .XS(tk)S t (t k ) 
k=l k=K 0 



(17) 



so oil : 

- B(t k ) est la contribution du bruit et des brouilleurs sur X(t k ), 

- S(t k ) est la contribution du signal utile sur Xfo), 

- Ko est I'instant d'apparition du signal utile sur le palier, 

- e est un terme prenant en compte rintercorreiation entre le signal B(t) et le signal utile S(t), qui tend vers 
55 0 lorsque K devient grand. 

Pour ia position de synchronisation testee : 



6 
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fXd = K 



K 0 +K 

2 x(tk)d-(k) 



(18) 



10 



15 



20 



25 



' Dans les conditions asymptotiques (avec K»1), lorsque le signal utile est synchronise avec la r6plique 
d(t), Rxx et f X d se simplif ient et s'Scrivent : 



: R+ K '"^ + 1 n,SS' (19) 



A 



j-V^S (20) 

Le critere obtenu en estimant la matrice de correlation Rxx sur K* echantillons poss6de comme valeur 
asymptotique, apres un calcul utilisant le lemme d'inversion matricielle: 



C' 0 = 



SNIRopt 



! + K'-Ko + 1 SN|Ropt 



(21) 



k: 



Comme K * * 1 < 1 , alors C > > C 0 ' 



La valeur du critere apparait dans ces conditions augment^ d ; un facteur dependant de I'instant d'appa- 
rition du signal utile sur la frequence consid6r6e, ce qui entrame une amelioration de la prise de synchronisa- 
tion, au prix d'une puissance de calcul reduite. 

□'autre part, le nombre K' ne doit pas etre choisi trop grand, sa valeur est d&erminee pour que sur 
K* echantillons, I'environnement de brouillage puisse etre considere comme stationnaire. 

Le seuil e prendre en consideration depend de la probability de fausse alarme (Pfa). La (Pfa) est la pro- 
babilite qu'un bruit soit detecte comme une synchronisation. Le calcul de la valeur du seuil r\ a lieu en calculant 
Pfa=P(C>r|) lorsque Tantenne recoit uniquement du bruit. 
En presence de bruit, le critere de synchronisation s'ecrit: 



30 



C = 



» £ B K V 



^ K 



\-1 



77 2 B k B k* 
r k=1 



^I 8 ^ (22) 
K k=1 J 



35 



En supposant que B k est un bruit gaussien complexe de puissance o 2 , independant entre les capteurs, 
lorsque.K est suffisamment grand, le terme du milieu de la relation (22) devient : 



40 



/ r 1 K 
77 2 B k B k 



= cx +2 



Id (23) 



45 



(Id etant la matrice identite ou unitaire) 

le critere de synchronisation peut alors etre approxime de la facon suivante ; 



(24) 




50 



Soit: 



55 



1 " 
K a i=i 



k=1 



2^ 



K 2 rr 2 
K ° i=1 
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avec: 



K 



Z B k^ = Ri + Jli (26) 
k=1 



Bi represents ie bruit obtenu sur le capteur i, pour I'echantillon k. 

R, + jl, peut 6tre apprpxime par un bruit gaussien complexe de variance K o 2 . Dans ces conditions R t l t 
*0 sont chacun des bruits gaussiens r6els de variance 



15 



En notantJi = Rj ou l (l C devient : 



(Ko 2 ) 



20 



Cs ^ll.(f) (27) 



25 



30 



K a 2 

ou Tj est une variable aleatoire qui suit une loi normale centree de variance — — • La variable aleatoire Z\ 

i T| suit done une loi normale centre de variance 1 . 
VKo2 

C devient : 



1 



2N 



C£g2f (28) 



35 



Comme la suite des Zj peut etre considers comme une suite de 2N variables al6atoires independantes 
gaussiennes de moyennes nulles et de variances 1, la variable aleatoire 



40 



2N 

i=1 



45 



suit une loi du x 2 a 2N degr6s de Iibert6. 
La probability d'avoir C ^ n s'6crit alors : 



50 



p(C>n) = p 



2N ~ 



£zf >2K n 

t>1 N 



= Pfa (29) 



Le tableau suivant donne differentes valeurs de seuils r| obtenues par recherche dans une table donnant 
P(Z^A) pour une loi du % 2 a 2N degres d_e liberte et pour Pfa = 10- 6 . 
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N = 5 


=> 


A = 42 


=5> 


rx = 0,26 


N = 4 


=> 


A =37 


=> 


7i = 0,23 


N = 3 


=> 


- A = 33 


=> 


t| = 0,21 


5 N = 2 


=> 


A = 28 


=> 


r> = 0,18 



Comme entre remission et la reception, le signal transmis peut etre affects d'un certain decalage en fre- 
quence, ou decalage Doppler, du essentiellement aux derives en frequences des oscillateurs locaux des 
10 synthetiseurs a remission et a la reception, ainsi qu'aux conditions de propagation ionospherique, il est ne- 
cessaire de compenser ce decalage. Le risque est s'il n'est pas compense, d'empecher toute prise de syn- 
chronisation. 

En appelant Af 0 . le decalage Doppler entre le signal recu par Pantenne et le signal attendu et en considerant 
la position de synchronisation, le vecteur signal recu peut s'ecrire avec les notations des paragraphes sui- 
ts vants : 

X(t) = d(t) e^o-t S + B(t) 
La matrice de correlation estimee peut aiors s'ecrire, : 



20 



RXX = ^ Z [V^ dik)* 2 * "*** S + B(tk)J^ d(t k )e 2 -j- Af o ^.S + B(t k )|. (30) 



25 



40 



45 



K 

RXX = e + 7 l[7c s |d(t k )| 2 S S* + B(t k ) B(i k ) { ] (31) 
K k=1 L J 



ou e est une matrice prenant en compte les interrelations entre la contribution du signal utile et le bruit b(t). 
30 Lorsque K»1, dans les conditions asymptotiques, Rxx peut s'approximer par : 

Rxx^sSS* * cr2 Id (32) 
La relation (32) montre que le decalage Doppler ne modifie pas le calcul de Rxx- En effet, le decalage 
s'elimine naturellement lorsque le signal utile est multiplie par sa quantite conjuguee. 
Le vecteur d'intercorrelation peut s'ecrire : 

35 



ta-^ijVS^V^^S + ^Jd-Cg (33) ( 



f Xd = + ± i[^|d(t k )f ***** S] (34) 



oil e' est upi terme prenant en compte Tintercorrelation entre le bruit B(t) et la replique. Lorsque K»1 , r M peut 
s'approximer par : 



50 



1 K - 2jn.^a.k 

* k=l 



S (35) 



55 La relation (34) montre que le decalage Doppler modifie I'estimation du vecteur d'intercorrelation r* d et 

que du fait que la phase du terme en exponential complexe peut tourner de plusieurs fois 360' sur Hntervalle 
de temps utilise pour le calcul de la correlation, Integration ne peut plus avoir lieu de facon coherente ce qui 
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10 



15 



emp§che le critere de synchronisation muiticapteurs de converger. Le processus suivant permet d'estimer le 
d6ca!age Doppler et de le compenser pour rendre possible la prise de synchronisation. 

Selon ce processus, les K s/mboles de la s6querce de synchronisation sont decoup6s de la fa$on repre- 
sentee d la figure 5 en M tranches de Ko symboles (M K 0 = K), de telle sorte que sur chacune des tranches 
de Kosymboles, le d6calage Doppler fasse tournerla phase du signal re$u parl'antenned'une valeur inferieure 
d 90°, af in de permettre une integration coh6rente de r Xd sur chacune des tranches. 

Si r X di d6signe la valeur de rintercorr&ation estim6e sur la tranche i, celle-ci s'6crit : 



r Xdi = 



(M)K 0 2jn Afe k 

i7 s fern?* 

k=iK 0 +1 



S i = 0 M-1 (36) 



En mettant en facteur dans la relation (36) un terme dependant uniquement de Ko, la longueur de la tranche 
6l6mentaire consid6r6e, et du d6calage Doppler, la relation (36) devient : 



20 



r Xdi = 



j^Z>/^K(k + iKo)T)| 2 e F « 



2j7t 



iK 0 



• S i = 0 M-1 (37) 



25 



30 



Ce qui peut encore s'ecrire : 

f x<(l s x, e«» ^«. .S i = 0.....M - 1 (38) 
La relation (38) montre que les vecteurs rxdi sont des vecteurs, echantillonnfe a Leur transformee de 
Fourier discrete est alors : 



35 



40 



45 



M-1 



-2j*< 



Af 



rxdUf ) = Z r X di e Fe/K °' 
i=0 



fxd(Af)= X>qe 
i=0 



M-1 2j7t.^ _Af ■■ 



Fe 



iKo 



(39) 



(40) 



La relation (40) montre que lorsque Af = Af 0 , le d^calage Doppler est compens6, ce qui rend maximum la 
valeur du critere de synchronisation r^AO^xx^rxd^O- 

Af in d'estimer Af, c'est-S-dire af in d'estimer le decalage Doppler qui conduit £ la valeur maximale du cri- 
tere, et af in de rSduire le cout de calcul correspondant, une transformee de Fourier rapide est effectu§e sur 
les vecteurs r^i (le nombre de tranches el^menta ires'M etant choisi pour former une puissance de 2) : ce calcul 
donne : ' 



50 



M-1 



M-1 



rxdM^X r Xdl e . £ f Xdi e 

i=0 j=0 



-2jn* 



k = 0 M-1 (41) 



oli Afi est la r6solution Doppler telle que : 

Af , = f^2 (42) 

Si la resolution est insuf f isante, une amelioration d'un facteur 2 peut Stre obtenue en calculant une trans- 
formee de Fourier rapide, sur les vecteurs r Xdi "decalees" en frequence de la fa?on suivante : 



10 
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20 



f Xdi e M 



e M k = 0 M-1 (43) 



L'op6ration consiste a multiplier chacun des vecteurs d'intercorrelation r Xdi par le terme e-U"" M ) et a effec- 
tuer ensuite un calcul de transformee rapide sur les vecteurs ainsi obtenus. 

Le critere de synchronisation est ensuite estime sur chacune des cases de calcul des transform6es de 
10 Fourier pr6c6dentes par la relation : 

i*f)^fy-^f) ^0 M-1 («4) 

Bien evidemment, si. la resolution est.toujqurs.insuffisante, il est possible de Tam6liorer en utilisant le 
meme processus pour calculer des vecteurs ; 



r x ^KAf 1 + ^ 



ou n est un entier quelconque. 
25 Un test sur la case de la transform^ correspondant au maximum du critere de synchronisation est effeo 

tu6. Si le seuil est depasse pour une valeur ko y donn6e du decalage Doppler, II est alors possible d'affiner 

le calcul du dtcalage Doppler (indispensable pour la demodulation) en calculant la valeur du entire pour dif- 

ferentes valeurs de decalage Doppler entourant ko ^ En effet, k Tissue de cette etape, la resolution sur Pes- 

30 " 2 

timation du decalage Doppler est Pour I'ameiiorer, il faut proceder par dichotomie : k chaque iteration et 

de la facon representee & la figure 6, la valeur du critere pour la valeur de decalage Doppler situee au milieu 
des deux valeurs de decalage Doppler donnant les valeurs maxima du critere lors de Iteration precedente 
35 est calcuiee. La recherche par dichotomie est arret6e des que la resolution souhaitee sur I'estimation du de- 
calage Doppler est obtenue. Apres n iterations de I'algorithme de recherche par dichotomie, la resolution, qui 

etait initialement de est : 

40 ~zT - 2 n + i < 45) 

. L'enchainernent des 6tapes de calcul pr6c6dente pour la mise en oeuvre du proc6d6 selon ['invention k 
I'aide de moyens de traitement appropries form6s par example par des processeurs de traitement du signal 
ou tout dispositif equivalent connu est f igur6 sur Porganigramme de la figure 7. 

II consiste suivant I'etape 9 £ estimer la matrice de correlation Rxx suivant la relation (7). Cette matrice 
45 est estimee sur un bloc contenant un nornbre d'echantillons K f superieur au nombre d'echantillons de la se- 
quence de, synchronisation. 

[.'estimation de Rxx est utilis6e ensuite pour toutes les positions de synchronisation incluses dans le bloc. 
Puis la matrice de correlation inverse R« 1 est cateuiee a I'etape 10. f Xdi qui estdefini par la relation (8) est cal- 
cul6 k I'etape 11. Pour chaque 6chantillon correspondant k une position de synchronisation P dans la suite 
50 des echantillons P = 1 k P = K'-K, f Xdj est calcuie pour i = 1 k M sur M blocs de Ko symboles comme d6f ini k 
la figure 5 suivant la relation 
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rxdi = s?- Z X(P+K + iK 0 )d (K + iKq) (46) 
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Des calculs de transform6e de Fourier Rapide sont effectu6s & l'6tape suivante 12 sur un nombre deter- 
mine M (M = 8 par exemple) de vecteurs f X( ji suivant la relation (41). 

Aux Stapes 13 et 14 un calcul de vecteur R Xd i d6ca!6 et un calcul de transformee de Fourier rapide suivant 
la relation (43) sont executes pour ameiiorer la resolution par un facteur de 2. 

Le crittre de synchronisation est calcuie k I'etape 15 en application de la relation (44). 

La comparaison du critere ou seuil r\ est effectuee a I'etape 16. Si a I'etape 16 la valeur du critere 

est superieure ou egale au seuil, la position de synchronisation est d6finie a I'etape 17 sur la position P de 
I'echantillon correspondant. L'estimation du d6calage Doppler peut ensuite etre aff in6 & I'etape 1 8 en calculant 
a nouveau la valeur du critere de synchronisation pour des d6calages Doppler entourant le d6calage 

ho 

precedent. Cette etape peut etre suivie par une etape 19 d'estlmation de la position de synchronisation pour 
les autres trajets. 

La recherche de synchronisation prec6demment decrite est effectu6e pour chacun des 6chantillons de 
signal num6ris6. En pratique, pour ameiiorer la precision de la position de la synchronisation il est preferable 
d'effectuer un sur6chantillonnage du signal d'un facteur determine L, L £tant 6gal a 4 par exemple. Les trai- 
tements de Torganigramme de la figure 7 et dScrits ci-dessus sont encore utilisables sur le signal ainsi sure- 
chantillonne, mais il est possible d'utiliser une variante dans le proc£d£, af in de r6duire la puissance de calcul 
n6cessaire au traitement. Dans ces conditions, aprfcs avoir calcuie pour chaque echantillon le vecteur r Xd i pour 
i « 1, ... M sur M blocs de Ko symboles suivant la relation (46) un entire de synchronisation C ( est calcuie sur 
chacun des M blocs par la relation 

C, = ^ xdr R^r Xdl (47) 
Puis un critere somme de I'ensemble des C, est calcuie suivant la relation : 

Ci(P) = £C, (48) 
i=1 

Ce critere correspond a une integration non coherente des Ci (contrairementa ce qui est fait pour le calcul 
par FFT ou on intfcgre de facon coherente), ce qui se traduit par une remontee du crit£re sur les echantillons 
de bruit. 

Le calcul du critere a lieu par blocs de L echantillons cons6cutifs. Sur chacun de ces blocs, I'echantillon 
qui donne la valeur maximum de C^P) est s6lectionne pour effectuer le calcul par transforms de Fourier ra- 
pide necessaire a I'estimation du decalage Doppler comme deja decrit aux Stapes 11S 19 de la figure 7. 



Revendic^tions 

1. Precede permettant a un modem de se synchroniser sur un transmetteurdesignaux numSriques parvoie 
hertzienne en presence de brouilleurs consistant a : 

- estimer (9) la matrice de correlation Rxx de signaux recus sur un ensemble de N capteurs de recep- 
tion (1, ... 1 N ) par blocs composes d'un nombre d'echantillons d'une sequence de synchronisation 
transmise par le transmetteur, A 

- & calculer (10) la matrice de correlation inverse R^ 1 , 

. a calculer (11) des vecteurs d'intercorreiation f Xd entre les signaux X(k) recus sur I'ensemble des N 
capteurs et un signal de r6plique d(k) connu, 

- a calculer un critere de synchronisation (7) effectuant le produit scalaire entre le vecteur d'intercorre- 
iation r Xd et le vecteur obtenu par produit de Tinverse de la matrice de correlation R^ et du vecteur 



12 



EP 0 665 665 A1 



d'intercorreiation r X d 

- a comparer (8) la valeur du critere obtenu & une valeur de seuil y\ d6termin6e pour placer la synchro- 
nisation sur l'6chantillon pour lequel la valeur du critere depasse la valeur de seuil r\ t caracterise en 
ce qu'il consiste pour estimer le d6calage Doppler inherent aux communications entre 6metteurs et 
recepteurs, a d6couper la sequence de synchronisation en un nombre determine de M tranches de 
k 0 symboles, a effectuer (11, 13) surces tranches des calculs de vecteurs d'intercorreiation partielles 

f Xdi (i = 1 M) de telle sorte que sur chacune des tranches de ko symboles, le decalage Doppler 

fasse tourner la phase d'une valeur inf6rieure a 90°, a effectuer (12, 14) des calculs de transform^ 
de Fourier rapide (FFT) sur les vecteurs f X di obtenus, et a calculer une valeur du critere de synchro- 
nisation a partir des vecteurs ^(kA^ obtenus en resultat des calculs de transform6e de Fourier ra- 
pide, et de la matrice de correlation Rxx- 

A 

Proc6d6 selon la revendication 1 caract6ris6 en ce que ('estimation de la matrice de correlation Rxx est 
obtenu par la relation ' 



ou p est la position du premier echantillon a partir duquel la correlation est effectu6e entre le vecteur X 
et la replique d. 

Proced6 selon la revendication 2 characterise en ce qu'il consiste a choisir la case de calcul de la trans- 
formee rapide qui donne une valeur maximum du critere de synchronisation et a comparer la valeur maxi- 
mum obtenue a une valeur de seuil determinee pour decider de la synchronisation. 

Proc6d6 selon Tune quelconque des revendications 1 a 3 caracterise en ce qu'il consiste pour obtenir une 
meilleure resolution sur le d6calage Doppler et une valeur plus importante du maximum du critere, a cal- 
culer (13) des vecteurs r xd , decal6s en frequence et a effectuer un calcul de FFT sur les vecteurs r xd , de- 
cays de facon a obtenir la valeur du critere avec une meilleure resolution sur les frequences Doppler. 

Precede selon Tune quelconque des revendications 1 a 4 caract6rise en ce qu'il consiste a sur6chantil- 
lonner les signaux repus sur les N capteurs par rapport au rythme symbole, pour ameiiorer la resolution 
sur la position de synchronisation. 

Procede selon I'une quelconque des revendications 1 a 5 caracterisee en ce qu'il consiste a effectuer le 
calcul de la matrice de correlation Rxx sur un nombre d'echantillons K* superieur au nombre K d'6chan- 
tillons de la sequence de synchronisation af in de diminuer la puissance de calcul. 

Proc6d6 selon Tune quelconque des revendications 1 a 6 caracterise en ce qu'il consiste pour diminuer 
la puissance de calcul et apres avoir sur6chantillonn6 d'un facteur L par rapport au rythme symbole, a 
effectuer une etape de pretraitement consistant a calculer pour chaque echantillon p (p = 1 a L) composant 
un symbole les valeurs C, du critere de synchronisation sur chacune des M tranches de k<, symboles et 
a calculer la valeur du critere en effectuant la somme des criteres C, obtenus pour retenir dans la suite 
du traitement rechantillon donnant la valeur maximale du critere C = ICj. 

Dispositif pour la mise en oeuvre du proc6d6 selon Tune quelconque des revendications 1 a 7 du type 
' comportant, un r6seau de N capteurs (1 ,...1 N ) de signaux de reception couple a un r6seau de f litres adap- 
t6s (2) comportant des circuits de pond6rations amplitude phase commandes par un processeur de 
commande (3), le processeur de commande (3) etant programme pour 




et le vecteur d'intercorreiation est estime suivant la relation 




K 
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estimer la matrice de correlation inverse R« des signaux repu X (K) et calculer les vecteurs d'inter- 
correlation r Xd entre les signaux X (K | repus etle signal de replique d (K) caracterise en ce qu'il comprend 
des moyens de traitement (5... 19) pour calculer un critere de synchronisation (7) en effectuant un 
produit scalaire entre le vecteur d'intercorr6lation r xd et le vecteur obtenu par le produit de I'inverse 
de matrice de correlation R« et du vecteur d'intercorr6lation ?xd. 

comparer (8) la valeur du critere obtenu & la valeur de seuil r\ pour placer la synchronisation sur 
i'6chantillon pour lequei la valeur du critere d&passe la valeur du seuil f), 

et estimer le decalage Doppler inherent aux communications entre 6metteurs et recepteurs, en d6- 
coupant la sequence de synchronisation en un nombre determine de M tranches de k<, symbole^, 
en effectuant (11, 13) sur ces tranches des calculs de vecteurs d'intercorr6lation partielles r xdj (i = 

1 M) de telle sorte que sur chacune des tranches de Kq symboles le decalage Doppler fasse tour- 

ner & la phase d'une valeur inf6rieure & 90°, en effectuant (12, 14) des calculs de transform6e de 
Fourier rapide (FFT) sur les vecteurs r Xd| obtenus, et en calculant une valeur du critere de synchro- 
nisation & partir des vecteurs r Xd (KAf]|, obtenus en r6sultat des calculs de transform6e de Fourier ra- 
pide et de la matrice de correlation Rxx* 
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